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Abstract  —  This  paper  is  concerned  with  a  quantitative  nondestructive  evaluation  of  conduc¬ 
tors  using  superconducting  quantum  interference  devices  (  SQUIDs  ).  A  measurement  system 
is  described  for  an  electrical  potential  problem  with  an  unknown  boundary.  Domain  identifi¬ 
cation  is  discussed  within  the  theoretical  framework  of  a  parameter  estimation  problem  for  the 
electrostatic  field  analysis.  Applying  the  method  of  mappings  to  the  problem  considered  here, 
we  present  computational  methods,  including  theoretical  convergence  results  for  the  associated 
finite  dimensional  problem  identification  techniques. 

1.  INTRODUCTION 

Recently,  demand  has  grown  for  assessing  the  structural  integrity  of  materials  used  in 
nuclear  power  plants  using  advanced  sensor  technologies.  An  important  effort  on  such 
problems  entails  quantitative  nondestructive  evaluation  methods  in  magnetic  flux  imag¬ 
ing  techniques.  These  methods  involve  an  attempt  to  characterize  structural  flaws  or 
defects  which  may  not  be  detectable  by  visual  inspection.  SQUIDs  (  superconducting 
quantum  interference  devices  )  have  the  potential  to  detect  material  defects  in  conductors 
due  to  their  extremely  high  magnetic  flux  sensitivity  [1],  In  this  paper,  we  propose  a 
computational  method  for  recovering  defect  shapes  with  magnetic  flux  density  data  from 
high  critical  temperature  SQUIDs  (HTc-S QUIDs). 

In  the  proposed  nondestructive  test,  a  stationary  current  density  is  applied  to  the  conduc¬ 
tor  inspected.  The  magnetic  flux  density  can  be  measured  from  a  SQUID  sensor  located 
near  the  conductor.  Figure  1  illustrates  the  inspection  process  using  an  HTc-SQUID.  As 
shown  in  Fig.  1,  the  existence  of  a  defect  (  represented  by  a  nonconducting  volume  )  cor¬ 
rupts  the  current  flows  inside  the  conductor  and,  as  a  result,  this  material  defect  can  be 
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detected  as  the  perturbation  of  magnetic  flux  data.  It  is  well-known  that  a  mathematical 
model  for  this  nondestructive  test  can  be  derived  from  Maxwell’s  equations.  By  Faraday’s 
law  and  Gauss’s  law,  the  electrical  potential  E  satisfies 

V  x  E  =  0  (1) 

V  ■  eE  =  p  (2) 

where  e  and  p  denote  the  permittivity  and  the  charge  density  of  the  sample  specimen, 
respectively.  We  introduce  the  electrical  scalar  potential  (p  such  that 

E  = 


Fig.  1  Inspection  process  using  HTc-SQUID 


Assuming  that  there  is  no  charge  density  inside  the  material  and  that  the  conductor  is  a 
homogeneous  plate,  we  can  rewrite  Eqs.  (l)-(2)  as  the  Laplace  equation 

VV  =  0  in  G3C  R3,  (3) 


with  non-homogeneous  boundary  conditions 


cf)  =  v  on  dG!/.  V<f>-h  =  0  on  OG*  U  !)(% . 


where  G 3  and  dG3  =  dGfudG^  UdG 3  denote,  respectively,  the  domain  and  the  boundary 
of  the  conductor  with  a  flaw.  More  specifically,  the  boundary  dGf  U  dG 3  denotes  the 
exterior  boundary  of  the  conductor  and  dG({  characterizes  the  defect  shape.  In  the  above 
equations,  v  denotes  the  prescribed  electrical  potentials  on  dGf .  The  measurement  of 
magnetic  flux  density  by  HTc-SQUIDs  is  described  by  Biot-Savart’s  law  [2],  for  x  = 
(•Ti,  x2,  x3)  G  R3  -  G3, 


B(x) 


cr 0P0  r  v<f)  X  (x  -  .-r')  , 

47T  Jg3  I  x  —  x'  |3 
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where  <r0  and  //0  denote  the  electrical  conductivity  and  magnetic  permeability  of  the 
conductor,  respectively.  Such  measurements  can  be  obtained  by  detecting  voltage  changes 
in  SQUID  magnetmetors.  We  assume  that  the  measurements  are  only  made  on  the  vertical 
direction  (.r3  axis).  Since  the  voltage  change  can  be  measured  through  the  pickup  coil, 
we  can  only  take  the  average  of  the  magnetic  flux  density  at  the  sensor  location.  Suppose 
that  ip  £  R3  —  G3  is  the  location  of  sensor  and  Sp  denotes  the  region  of  the  pickup 
coil  with  the  center  xp.  The  problem  considered  here  is  that  of  identifying,  from  input 
and  output  data  {v.  /i(.r  ..)},  the  geometrical  shape  of  the  defect.  If  we  assume  that  the 
conductor  is  a  thin  plate  and  that  the  defect  to  be  identified  is  uniform  in  the  x3  direction, 
the  corresponding  electrical  potential  becomes 


Vcf)  ~ 


d(f> 

d^2 


(4) 


Then  the  observation  at  xp  can  be  approximated  by 


BnM  *  v,  =  |/0j  {(*»  -  4)§*r  ~  (*•  -  *i)|| 


,,  d<f> 


\x  —  x 


j  1-3 


dx 


dx . 


(6) 

where  | |  denotes  the  crosssectional  area  of  the  HTc-SQUID.  From  our  assumptions,  the 
identification  problem  is  well  approximated  by  using  a  2-D  spatial  domain.  Thus,  in  this 
paper,  our  attention  is  restricted  to  a  domain  identification  problem  in  two  dimensions. 


This  paper  is  organized  as  follows.  In  Section  2,  the  mathematical  model  for  the  system 
is  described  by  the  Laplace  equation  in  a  two-dimensional  spatial  domain.  The  measure¬ 
ments  are  derived  from  Biot-Savart’s  law.  Then  these  problems  are  treated  as  domain 
identification  problems  in  electrostatic  held  analysis.  In  Section  2,  we  also  formulate  this 
problem  in  an  abstract  setting  in  a  Hilbert  space.  The  ideas  proposed  in  Banks  et.  ah, 
[3]  [4]  using  the  “method  of  mappings”  are  adopted  to  the  problem  considered  here.  In 
Section  3,  for  computational  purposes,  we  approximate  the  Hilbert  space  by  finite  dimen¬ 
sional  subspaces  and  we  discuss  convergence  analysis  for  the  approximate  identification 
problems.  A  practical  computational  algorithm  is  proposed  briefly  in  the  last  section. 


2.  PROBLEM  FORMULATION  AND  BASIC  ASSUMPTIONS 


Let  G  be  the  sectional  plane  of  the  sample  specimen  G3  as  shown  in  Fig.  1.  Let  q  be 
a  constant  vector  which  characterizes  an  unknown  defect  domain  where  q  belongs  to  an 
admissible  parameter  set  Q  such  that 

(H-0)  The  admissible  parameter  set  Q  is  a  compact  subset  of  RM . 

We  consider  a  bounded  domain  Cq ( C  G  C  R 2)  which  is  parametrized  by  unknown  values 
of  q.  As  depicted  in  Fig.  2,  the  bounded  domain  Gq  is  considered  as  the  system  domain 
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such  that  Gq  =  G  —  Cq.  Corresponding  to  the  boundary  0G!d .  dG^  and  0G(A  in  three 
dimensions,  we  assume  the  boundary  of  Gq  consists  of  the  following  components: 

dGq  =  u  UidGi  U  dCq 

We  further  assume  that  meas(Gg)  >  8  >  0,  uniformly  for  q  £  Q,  for  some  positive  8. 


Fig.  2  Spatial  domain  G  in  two  dimensions 


Thus,  from  Eqs.(3)  and  (4),  the  system  can  be  rewritten  by 

V2cf>  =  0  in  Gq  (6) 

(j)  =  vt  on  dGi,  i  =  1,3,  (7) 

V <f>  ■  h  =  0  on  dGi  U  dCq,  i  —  2,4.  (8) 


Suppose  that  the  SQUID  sensor  is  scanned  on  the  x1  —  x2  plane  and  measurements  are 
made  at  { x *p  =  (.r1,  ,  .r,^  )}"  Let  S‘p  be  the  region  of  the  pick-up  coil  whose  center  is 
located  at  xp.  We  also  assume  that  the  distance  between  the  surface  of  the  plate  and 
SQUID  sensor  is  taken  as  h.  Then,  from  (5),  the  observation  in  case  of  two-dimensions  is 
described  by 


^  = 


cr0fi0A  f  f  \  (J(p  Ocp 

\k  r*  (ar*  "**'53 


d(f> 


d(f> 


x  {(xj  -  x'Q2  +  (x2  -  x'2)2  +  h2} 
for  i  —  1,2, 


-3/2 


d.Xq  dx'2 


dxidx  2, 


,  m . 


(9) 
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where  A  denotes  the  thickness  of  the  plate.  Let  0*  be  an  arbitrary  smooth  function  on 
Gq  satisfying 

(p*  =  Vi  in  dGi ,  i  =  1,3. 

Then  the  system  (6)- (8)  is  equivalent  to  the  following  boundary  value  problem  : 

0  =  U  +  0* 

where  u  is  the  solution  of 

V2u  =  -VV  in  Gq  (10) 

u  =  0  on  dGi,  i  —  1,3,  (11) 

Vv/  •  n  =  'V&*  ■  ti  on  ()G.j  U  ()C,..  i  —  2,4.  (12) 

Let  Vq  be  the  closed  subset  of  functions  in  H1(Gq)  satisfying  the  homogeneous  boundary 
conditions  on  dGi ,  *  =  1,3,  i.e., 

Vg  —  {  0  |  0  £  H1(Gq),  0  =  0  on  z  =  1,3  }. 

endowed  with  the  norm 

mVq=  [  m\2dx,  (is) 

JGq 

Lemma  1.  Suppose  that 

(c-i)  e*  is  an  element  of  V'\  the  inhomogeneous  analogue  of  Vq  dehned  by 
Vq  =  {  0  |  V  €  // 1  (G,.j,  0  =  Vi  on  <9G';,  *  =  1,3  }. 


Then,  for  every  q  (E  Q,  there  exists  a  unique  solution  (f>(q)  G  V*  of  (6)-(8)  in  the  sense 
that 

U  =  (j)  —  <f>*  <E  Vq 

satishes  ( 10)- ( 12).  The  solution  u  of  ( 1 0)-(  12)  is  in  the  weak  sense,  i.  e., 

a(q)(v,  ih  •  a{q)(&*,  v)  for  0  G  Vq.  (14) 

Here  a(q){-.  •)  denotes  a  sesquilinear  form  on  Vq  X  Vq  dehned  by 

cr(q)((p,  0)  =f  f  Vip-'Vipdx.  (15) 


Moreover,  we  have 

N<z)k  <  Ki  W\Vq 

where  K-\  is  a  constant  independent  of  q. 
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Remark.  We  note  that  even  though  (if  is  in  V*  and  not  in  Vg,  the  value  of  a(q)(cf.  tjf) 
for  0  £  Vq  can  be  defined  as  in  (15)  and  the  value  of  is  well  clehned  as  in  (13).  We 

tacitly  assume  in  the  statement  of  Lemma  1  and  in  subsequent  discussions  that  such  an 
interpretation  is  understood. 


Proof.  From  (15),  for  arbitrary  <p  £  Vq .  cr(q)  is  V-elliptic  with  constant  C\  =  1  which 
is  independent  of  q  since 

=  t  |Vy>|  2dx  =  |V^|^2(Gg)  =  \(p\lq.  (16) 

JGq 

Since  we  have 

k('/)(rW)l  <  W\vq\^\vq- 
the  linear  functional  a(q)(df.  •)  is  bounded  on  Vq ,  i.  e., 

I  -  I  <  for  9?  £  Vq.  (17) 

From  (16)  and  (17),  we  can  apply  the  Lax-Milgram  lemma,  taking  the  Hilbert  space  Vq, 
sesquilinear  form  a(q)(-,-),  and  linear  functional  —  a(q)((f>*:  ■),  respectively.  Hence,  for 
each  q  £  Q,  there  exists  a  unique  solution  u  £  Vq  in  the  sense  of  (14).  Similarly,  from 
(14),  (16)  and  (17),  we  have 

Mve  ^  I  “  cr(?)(^w)l  <  W\vq\u\vq 

from  which  follows 

\u\vg  < 

thus  completing  the  proof. 


For  convenience  of  theoretical  developments,  we  use  the  polar  coordinate  system  x  =  (r,  0) 
instead  of  the  Cartesian  coordinate  x  =  (x1,x2)  in  the  sequel.  To  this  end,  we  introduce 
the  unknown  defect  function  c(6,q)  and  the  prescribed  outer  boundary  function  / f  0 ) . 
These  functions  are  assumed  to  be  defined  on  (0,  2tt)  and  are  27r-periodic,  respectively. 
As  depicted  in  Fig.  3,  it  is  further  assumed  that  the  system  domain  Gq  can  be  described 

by 

G9  =  {  (r,  6*)  |  0  <  6*  <  27 r,  c{6 ,  q)  <  r  <  1(6 )  }, 

and  hence  6  — >  c(6,  q )  is  a  parameterized  function  which  is  assumed  to  characterize  the 
unknown  defect  shape.  The  boundary  of  Gq  is  also  defined  by 

dGi  =  {  (r,  6)  |  6i_  i  <6  <  6i ,  r  =  1(6)  }  for  i  =  1,2,  3,4, 
dCq  =  {  (r,  0)  |  0  <  6  <  27T,  r  =  c(6 ,  q)  }, 


where  60  =  0  and  04  —  2tt.  respectively.  The  sesquilinear  form  (15)  can  be  rewritten  as 

r2ir  fl(6)  (  dip  dip  1  d(f  d(Jj\ 


nlyo) 

( 6,q ) 


r  dr  dr  +  r  d6  d6  )  ^ ' 


(18) 
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Fig.  3  The  spatial  domain  G  and  its  boundary  using  the  polar  coordinate  system 

For  the  discussions  here,  we  restrict  the  geometrical  structure  of  the  boundary  dGq  by 
imposing  the  following  hypotheses: 

(H-l)  For  each  q  (z  Q,  c(-,q)  £  14^(0,  2tt)  and  1  £  l/F^o(0,27r). 

(H-2)  There  are  constants  H  and  l32  such  that,  for  q  (E  Q, 

0  <  /h  <  c(0.  q)  <  l32  <  1(6)  <  oo  for  6  in  (0,  27r). 

(H-3)  There  exists  a  function  d,  :  Q  X  Q  — >  R1  with  d,(q.  q)  — >  0  as  \q  —  <?|  — >  0  such  that 

|c(',?)  -  c(-,g)|i,oo  <  d(q,q)  for  q.  q  G  Q, 
where  |  •  |ljQO  denotes  the  norm  ofwy. 

Following  standard  procedures  in  the  method  of  mapping  techniques  ([3], [4], [9]),  we  in¬ 
troduce  the  affine  mapping  x  =  T(q)x  where  x  =  (r,  9)  and  x  =  ( r,6 ).  The  precise  form 
of  T(q)  in  this  paper  is  given  by 

|  r  =  j(q]f,6)  =  {c(Q,q)  -  l(6)}{c- -  1(6)}  +  1(6) 

I  6  =  6 

Let  G  be  the  reference  domain  given  by 

G  =  {  x  =  (r ,  6)  |  0  <  6  <  27t,  c  <  r  <  1(6)  } 
which  is  independent  of  the  parameter  q  as  depicted  in  Fig.  4. 
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W) 


Fig.  4  The  reference  domain  G 


Then,  under  this  coordinate  change  T(q)'j  the  reference  domain  G  is  transformed  into  the 
unknown  domain  Gq.  Let  us  define  the  Hilbert  space  defined  on  G, 

V  =  {  tp  |  0  £  fh1(G),  1} ')  =  0  on  dGi,  i  =1,3} 


endowed  with  the  norm 


Noting  that 


we  have 


2-7T  /*/($) 


dtp 

dr 


dtp 

~dQ 


drdd. 


VT(q) 

|VT(g)| 


d'jjdr  d'j/dd 
0  1 
d'y 

df'- ’ 


dp 

dr 


dp 

W) 

drdd 


£hy\  1  d(p 
dr  J  dr 

/chA  _1  dft 

[df)  ~d0^Jr~d0 


^ drdB . 
dr 


(20) 
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Hence  (18)  can  be  transformed  into  the  sesquilinear  form  on  V  X  V  as  follows: 


/  w  n  def  /27r  [lW  j  ,  .d(pdi/> 

=  l  l  W«)  Sf  gf  +  “’(«) 


dfS 


{  dip  dij)  dip  dil’ 
\dr  d9  +  d9  dr 


where 


«i(<z) 

02(g) 

03(g) 


-1 


X  +  1%!) 

1  d9  v 


=  -Mr%wr 


(21) 


Lemma  2.  With  the  hypotheses  (H-0)  to  (H-3).  there  exist  positive  constants 
a,  A,A2,  and  K3  such  that,  for  q,q  £  (J.  the  sesquilinear  form  d(g)(-,  •)  satisfies  the 
following  inequalities  for  all  ip,  xjj  £  V: 


where 


g  )(•</', '/•) 

> 

ly 

(22) 

^(gX^VOI 

< 

(23) 

l<%)(<^X)  -  V5)  1 

< 

(24) 

d(q.  q)  —>  0  as  \q  —  q\  — »  0. 


Furthermore,  a ,  Ii2,  and  /\3  can  be  chosen  as  constants  which  are  independent  of  the 
parameter  vector  q. 


Proof.  From  (21),  the  associated  quadratic  form  is 


ai(g)^i  +  2o2(g)66  +  «3(g  )C2 


^i2 


2dfde^7  + 


For  any  quadratic  form  QF  =  —  2 6^ £2  +  cC2  with  a  >  0,  c  >  0,  completion  of  the 

square  arguments  yield  QF  >  s2(ac  —  b2)/c  and  QF  >  ^2(ac  —  62)/g  so  that 


2 QF  > 


ac  —  b 2 
a 


t 


2 

2  • 
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Choosing 


2  ( ®l\2  i  di  ,  ( dl\2 

a  =  72  +  — ^  ,  b  =  r  and  c  =  — 

\d0  dr  d6  \dr 


we  easily  find  that 


ai  {clKi  +  2a2(?)^lS2  +  «3(<z)s2 


1  W 

-  27  Uf 


-1 


Note  that 


4l  +  laf I  I  +  va«. 


dj_  =  c(6,g)  -  1(6) 
dr  c  —  I  (O') 


-l 


3 


and  from  (H-2)  and  (19),  there  exists  a  constant  C2  independent  of  q  such  that 

di~*  A{c-m} 


7  dr  ~  fa- 1(0 ) 


>  c2 


where 


r  ■  li  {c-m) 

(_2  =  mm 


ee[ 0,2tt]  fij  —  1(0) 
With  (H-l),  there  exists  a  constant  C3  such  that 

S’MSIT- 


where 


c3 


min 

ee[0,2?r] 


82- 1(0)}  V 
c-m  } 


max  \l(0)}2  +  sup  sup 
oe[o,2i<}  geQ  s-ed 


-1 


can  be  chosen  independent  of  q  (  here 


we  use 


(H-0)  and  (H-3)  ).  Thus  if  we  choose 


«  =  —  C2  min(l,  C3), 


we  obtain  the  coercivity  (22)  of  the  sesquilinear  form  cr(q)( •,  ■).  To  prove  the  boundedness, 
we  note  that 

\8W)(y-  W'l  <  4  sup  _  |(  • 

i<3,x£G 

From  (H-0)  -  (H-3).  we  have 


sup  \ai(q,x)\  <  K2  for  i  =  1,2,3 

zeG 
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where  I{2  is  some  constant  independent  of  q.  This  implies  the  boundedness  (23).  To 
establish  the  continuity  property,  we  note  that,  for  any  q,q  £  Q, 


\(J 


WA'P-.i')  - 


< 


(2*  (m\s  (\  n 


+Wl)-^)}W: 

1  n  36  36 


drd6 


Under  the  hypotheses  (H-l)  and  (H-3),  we  can  argue  that  «,;('/),  i  =  1,2,3,  are  continu¬ 
ous  in  Lco{G).  We  can  thus  infer  the  continuity  of  the  sesquilinear  form  (18)  with  respect 
to  q  £  Q  and  the  proof  of  Lemma  2  is  thus  completed. 


Lemma  3.  Let 

u  =  u  o  T(q),  cj>:  ■—  (f>  o  T(q),  and  3*  =  <f>*  o  T(q ) 

be  the  transformed  system  state,  the  transformed  solution,  and  the  transformed  function 
associated  with  </>*,  respectively.  Then,  for  every  q  £  Q,  there  exists  a  unique  solution 

(25) 

in  the  sense  that 

a(q)(u,i/>)  =  —b (<?)(</>*,  VO  for  y  £  V.  (26) 

Moreover  the  solution  u  in  the  system  on  G  is  bounded  in  V  uniformly  in  q  £  Q. 


Proof.  Using  Lemma  2  and  from  (18),  (21),  and  (20),  we  have  that,  for  ip  £  V„  and 
<P  =  V  o  T(q),t 


Similarly,  by  virtue  of  Lemma  2,  it  can  be  easily  shown  that 


M  2Vq<k2\p\l 

where  k2  does  not  depend  on  q  £  Q.  This  implies  the  V^-norm  is  equivalent  to  the  norm 
in  V  uniformly  in  q  £  Q.  Thus  we  have  the  solution  u(q)  =  u(q)  o  T(q )  in  the  sense  of 
(26)  satisfying 


IM 


(<z)|y<a  2W{q)\vq  <  a  2  <  a  2  /\'|L 


\v 
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which  follows  from  Lemma  1  and  verifies  the  claim  of  Lemma  3. 


Lemma  4.  With  the  hypotheses  (H-0)  to  (H-3),  q  — >  u(q)  is  continuous  from  Q  to 
V. 


Proof.  Let  qk  — >  q  in  Q  and  let  ii(qk),u(q)  be  corresponding  solution  of  (26).  That 
is, 

a(qk){u{qk)^)  =  for  ip  £  V  (27) 

a{q)(v(q).  </’)  =  -d(q)(^*,^)  for  f  G  h  (28) 

Substracting  (28)  from  (27),  we  have 

a(qk)(u(qk% $)  -  a{q)(u{q),  ^)  =  [<%)  -  a(qk)](4>*, #)  for  all  $  G  V. 

Then  we  may  write 

a(qk)(u(qk)  -  u(q ),  f)  +  |d(7/fcj  -  d(q)](u(q),  V>)  =  [<%)  -  d(qk)](f,  tj?)  for  G  V.  (29) 

We  note  that  while  <f>*  (fc  V  —  because  of  its  boundary  values  it  is  only  in  H1(G )  — 
it  is  readily  seen  from  the  proof  of  Lemma  2  that  (24)  holds  for  p>  =  <f>* .  Choosing  as 
yt  =  ii(qk)  —  u(q )  in  V  and  using  (22)  and  (24)  in  Lemma  2,  we  find  from  (29) 

a\ii(qk)  -  <>(</) |p  <  d(qk:q){\^*\y  +  \u(q)\y}\u(qk)  -  %)|^:} 

or 

a\ii(qk)  -  ii(q)\y  <  d(qk,q){\(f)*\y  +  \ii(q)\y}  - +  0  as  qk  q. 

This  yields  the  desired  continuity. 


From  (9),  the  system  output  has  the  form 

oo/ioA 


K)  =  - 


f  r  r*  r  {*.(«-«')  A -icos(« +  6')# 
Js>pJo  Jc(e,q)  1  1  ’dr'  r>  1  ’  86' 


47t|5||  Js’p  Jo  Jc(8,q)  \  dr 

X  | r2  +  r'~  —  2 rr'  cos (6  -f  O')  +  h2|  2  drdO  dr' dO' 

for  i  =  1,  2,  •  •  •  m.  (30) 

From  (30),  the  observation  output  for  the  transformed  system  can  be  written  as 

Y(q)  =  'H{q)&(q).  (31) 

where  'H(q)  :  V  — >  Rm  is  given  by 

[H{q)q?\i  =  -  4 m jsi  j  J  j/»i(r,0,f?3^';g)|^  -  h2(r,  6,  r,  & ]  g)|£j  drdOdr  dO' 


P I  P 
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where 


for  i  =  1,  2,  •  •  •  ,  ra, 


h1(r,9,r  ,9';q) 

h2{r,d,r,d'-,  q) 
h(r,9,r  ,e'-,q) 


h(r,9,r  ,6  \  q)[sin(0  —  9  ) 

+{l{r,  h  g)}"1^ ^  —  cos(0  +  fl')]f 

o9 

h(r,  9 ,  r',  f)':  g){7(r',  fl:  g)}-1^^  I gl  cos(g  +  fl'), 
r  [r2  +  {7(1  ,  9]  g)}2  —  2r^{f  ,  <7)  cos(0  +  0  )  +  h‘ 


3 

2 


Lemma  4.  With  the  hypotheses  (H-0)  to  (H-3),  the  mapping  q  — >  Y(q )  is  continuous 
from  Q  to  it'" . 

Proof.  Taking  into  account  that  the  operator  'H(q)  belongs  to  C(V]  l\  T )  for  each 
q  £  Q,  we  see  that  the  above  statement  follows  directly  from  Lemma  3. 

Thus  the  identification  problem  of  interest  is  defined  by  the  Lt-to-data  functional 

E(q)  =  \\H(q)i(q)  ~  Ur,  (32) 

where  { Yrj  j  denotes  the  given  observed  data.  The  identification  problem  can  then  be 
stated  as  follows: 

(IDP):  Find  q*  £  Q  which  minimizes  £(g)  given  by  (32)  subject  to  (25),  (26),  and 
(31). 

Theorem  1.  Under  the  hypotheses  (H-0)-  (H-4),  the  problem  (IDP)  has  at  least 
one  solution  q*  £  Q. 

Proof.  With  Lemmas  3  to  5,  the  above  statements  follows  from  the  compactness  of 

Q- 

3.  APPROXIMATE  IDENTIFICATION  PROBLEMS 

The  computational  scheme  we  propose  below  is  based  on  the  use  of  a  finite  element 
Galerkin  approach  to  construct  a  sequence  of  finite  dimensional  approximating  identifica¬ 
tion  problems.  To  approximate  problem  (IDP),  we  choose  a  sequence  of  finite  dimensional 
subspaces  HN  C  V  such  that 

| PNif>  —  ip\y  — +  0  as  N  — »  00  for  £  V  (33) 
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where  PA  is  the  orthogonal  projection  of  H  =  L2(G )  onto  II x  .  The  approximating 
system  is  defined  for  iiN  ( q )  £  H A  by 

^(<?)(hJV(<?),^JV)  =  for  ^’A  £  iP\  (34) 

From  (31),  the  observation  output  for  the  approximate  system  (34)  can  be  represented  as 

YN(q)=HN(q)iF(q)i  (35) 

where  TiN (q)  :  II N  — >  Pm  and  where  <frN(q)  =  hA  +  </>*.  Then  we  seek  to  solve  the 
approximate  identification  problem: 

(AIDP)*  :  Find  qN  £  Q  which  minimizes 

EN(q)=l-\HN(q)4‘N(q)-Ydf 
subject  to  the  finite  dimensional  system  (34),  (35). 

Lemma  5.  Let  qN  —+  q  £  Q.  Then  uN(qN)  — +  u(q)  £  V. 

Proof.  We  have 

a(qN)(uN(qN)^N)  =  -c^)(^*)  for  ^  £  HN  (36) 

a{q)(u(q),ij))  =  a(q)[<j>m,i!.')  for  f  £  V.  (37) 

Then  it  follows  that 

I uN(qN)  -  T/(v)lv  <  K<A(<hV)  -  PNu(q)  \v  +  I PNu(q)  ~  %)|f 
Hence  from  (33)  it  suffices  to  prove 

|bA  (qN)  —  PNu(q)\y - ■>  0  for  qN  — >  q  £  Q  as  iV  — >  oo. 

Taking  y':  =  in  (37)  and  subtracting  this  from  (36),  we  obtain 

d(^)(h«(^)5  ,j,N)  ~  a(q)(n(q),  * N )  =  [< J(q )  -  a(qN)M,  4'N)  for  ^  <E  HN. 
Furthermore  we  have 

a(qN)(uN(qN)  -  PNu(q),  f N )  +  d(<?JV)(PJVh(g)  -  %),  0*) 

+  [*(9*)  -  ^(9)](fi(4)^)  =  [*(9)  -  d(^)](^,  ^). 

Choosing  0A  =  AA  =  uN(qN )  —  PNu(q ),  we  hnd  that 

*(9*)(A*  Aw)  +  a(g*)(P"fi(,)  -  %),  Aw) 

+  [%*)  -  <%)](%),  A*)  =  [<%)  -  cr!v  v)]!oA  Av). 
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From  this  we  obtain 

a\AN\l  <  K2\PNu(q )  -  h(g)|t|AW|t  +  +  d(q,qN)\^\y\AN^ 

Consequently,  we  have 

a|AA  \y  <  I<2\PNu(q )  -  v{q)\y  +  d(qN ,  q)(\u(q)\y  +  \fr\y). 

Thus,  given  any  qN  — >  q  G  Q,  it  follows  from  (33)  that  AA  — »  0  as  N  — >  oo. 

Since  it  can  readily  be  shown  that  the  approximate  solution  oA  depends  continuously  on 
q,  solutions  exist  to  the  problem  (AIDP)^  for  each  N.  Our  convergence  results  for  the 
approximating  identification  problem  (AIDP)a  are  summarized  in  the  following  theorem 
which  follows  from  standard  arguments  (  e.g.,  see  [7]  ). 

Theorem  2.  Suppose  that  hypotheses  (H-0)  to  (H-3)  hold  and  let  q A  be  a  solution 
of  the  problem  (AIDP)a.  Then  the  sequence  { q A}  admits  a  convergent  subsequence 
{qNk}  with  qNk  — >  q*  as  N &  — >  oo.  Moreover,  q*  is  a  solution  of  the  problem  (IDP). 

4.  COMPUTATIONAL  METHOD 

It  is  rather  straightforward  to  develop  computational  methods  based  on  the  preced¬ 
ing  foundations.  For  example,  we  may  consider  a  linear  spline  approximation  of  a 
parametrized  function  for  c(q,6).  Let  SM  be  the  set  of  piecewise  linear  splines  (  see  [5] 
for  more  details  )  with  the  knot  sequence  =  {«27r/JVf  and  basis  elements  { /LW  j  wn. 
Then  we  approximate  the  unknown  defect  shape  function  by 

M 

c(q ,  0)  =  ^  q™B™{6)  for  6  G  [0,  2tt].  (38) 

i=0 

In  order  to  ensure  (H-0),  we  impose  constraints  on  Q  expressed  by 

Q  =  {  q  €  RM  \  Pi<qi<fai  =  li%--,  M  }  . 

It  is  clear  that  the  defect  shape  function  defined  by  (38)  satisfies  the  hypotheses  (H-2) 
and  (H-3).  Let  us  choose  U^_1{^>Ar})Y_1  as  a  set  of  basis  functions  in  V.  That  is,  for  each 
JV,  are  linearly  independent  and  ^NSpan{lpP}fL1  is  dense  in  //1  (G).  Then  the 

approximation  subspaces  can  be  chosen  as  H A  =  span{ij)f  }f=1.  An  approximate  solution 
can  be  defined  by 

N 

iN  =  aN-+r  =  z,w?i,?  +  y  (39) 

8  =  1 

where  the  coefficient  vector  wA  =  {wf  }  is  chosen  such  that,  for  j  =  1,2 ,  •  •  ■  ,  A, 

N 
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Hence  the  system  can  be  approximated  by  solving  the  linear  system 


LN(q)wN  =  fN(q ) 

where 

[LN(q) ];a  =  %)«,  $*)  for  i,  j  =  1,2 ,  ■  ■  ■  ,  JV, 
LfW (9)]i  =  -*(«)(**,  V’f)  ^  t  =  1, 2,  ■  ■  ■  ,  iV. 

The  corresponding  output  can  be  computed  as 

YN(q)  =  H”(q)w”(q)+g(q) 


(40) 


(41) 


where 

[HN  =f  [W(g)]i#f  for  *  =  1,  2,  ■  ■  ■  ,  m;  j  =  1,  2,  ■  •  ■ ,  TV, 

and  where 

=f  t?^)]#  for  *  =  1,  2,  •  •  ■  ,  m, 

respectively.  Associated  with  the  problem  (AIDP)^,  the  problem  is  to  find  the  parameter 
(E  Q  which  minimizes 

=  \\YN(q)  -  Y*f  (42) 

subject  to  the  finite  dimensional  system  (39),  (40),  and  (41).  Solutions  to  the  problem 
(AIDP)^  exist  for  each  N .  In  order  to  implement  the  identification  scheme,  it  is  necessary 
to  evaluate  the  gradient  of  the  cost  functional  (42).  Let  q N  be  an  optimal  solution  of  the 
problem  for  (42).  Then  a  necessary  condition  for  q N  to  be  optimal  is  characterized  by 

VqEN(qN)  -(q-qN)>  0  for  all  q  G  Q, 


where  Vg  denotes  the  gradient  of  EN(q)  with  respect  to  q.  From  (42),  the  gradient  of 
E^  (q)  becomes,  for  i  —  1,  2,  •  •  •  ,  M , 


[VqEN(qN)]< 


/  dKN(g) 

\  % 


w 


N 


+ 


dfN 
%  ’ 


dHN(<i),..N  ,  d9 

o  ^  "T"  Q  5 

Oqt  Oq% 


YN(q) 


(43) 


where  vN  is  the  solution  of  the  adjoint  system, 


[KN(q)fvN  =  [' HN(q)]T(YN(q )  -  Yd). 


With  (43),  numerical  optimization  techniques  are  readily  applicable  to  the  problem  for 
(42)  (  see  [6]  and  the  references  therein  ). 

5.  CONCLUDING  REMARKS 
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The  results  given  here,  while  applied  to  a  particular  domain  described  by  the  interior 
boundary  c(0.  q )  and  exterior  boundary  1(0).  are  rather  generic  in  nature  for  elliptic  prob¬ 
lems  with  unknown  domains.  The  techniques  used  here  (  and  the  necessary  assumptions 
on  the  problems  )  represent  a  rather  nontrivial  improvement  of  the  results  and  ideas 
presented  in  Chapter  VI  of  Banks/Kunisch  [7]  for  parameter  estimation  convergence  in 
elliptic  problems.  In  particular,  the  presentation  for  the  elliptic  system  in  [7]  requires  a 
nonvanishing  positive  coefficients  in  the  lowest  order  term  (  i.e. ,  c(a;)  >  c  >  0,  see  (1.3)  of 
[7])  in  L(q)u  =  V  ■  (aV«)  +  b  ■  'Vu  +  cu  that  can  be  replaced  with  the  techniques  used  here 
(  The  results  of  [7]  cannot  be  applied  directly  to  the  problem  we  treated  here).  While  the 
ideas  are  in  the  sprit  of  [7],  the  detailed  arguments  are  more  similar  to  an  elliptic  version 
of  the  hyperbolic  system  presentation  of  Chapter  5  of  [8]  (  compare,  in  particular,  (A3N) 
and  (H7)  of  pages  124  in  [8]  with  our  conditions  (33)  and  (24). 

The  generic  nature  of  the  results  presented  here  is  primarily  manifested  in  the  nature  of  the 
transformation  7  (4/)  :  G  — »  Gq  between  the  reference  domain  G  and  the  unknown  domain 
Gq.  As  long  as  this  map  satisfies  requirements  on  |VT(</)|  and  [VT(</)]_1  needed  above 
and  yields  a  transformed  system  with  appropriately  smooth  coefficients  cq(</)(z  =  1,2,3), 
the  arguments  remain  valid  even  if  the  unknown  domain  is  not  “annular”  in  nature  with 
unknown  boundary.  For  example,  as  shown  in  Fig.  5,  a  parallelpiped  (  slit  or  crack  line  ) 
shaped  damage  and  structure  should  be  amenable  to  the  arguments  given  here  (  in  this 
case  we  would  use  rectangular  coordinates  in  the  place  of  polar  coordinates  )  where  now 
q  parametrizes  the  slit  or  crack-like  damage  (  similar  to  the  considerations  in  [3]  [4]  on 
thermal  nondestructive  evaluation  ). 
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Gq 


Fig.  5  Material  with  a  defect  of  slit  type  and  its  reference  domain  shape 
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